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ABSTRACT 

Context. Until a few years ago, the amplitude variation in the photometric data had been limitedly explored mainly because of time 
resolution and photometric sensitivity limitations. This investigation is now possible thanks to the Kepler and CoRoT databases which 
provided a unique set of data for studying of the nature of stellar variability cycles. 

Aims. The present study characterizes the amplitude variation in a sample of main-sequence stars with light curves collected using 
CoRoT exo-field CCDs. 

Methods. We analyze potential stellar activity cycles by studying the variability amplitude over small boxes. The cycle periods and 
amplitudes were computed based on the Lomb-Scargle periodogram, harmonic fits, and visual inspection. As a first application of our 
approach we have considered the photometric data for 16 CoRoT FGK main sequence stars, revisited during the IRaOl, LRaOl and 
LRa06 CoRoT runs. 

Results. The 16 CoRoT stars appear to follow the empirical relations between activity cycle periods (Pcyc) and the rotation period 
(Prot) found by previous works. In addition to the so-called A (active) and I (inactive) sequences previously identified, there is a 
possible third sequence, here named S (short-cycles) sequence. However, recovery fractions estimated from simulations suggest that 
only a half of our sample has confident cycle measurements. Therefore, more study is needed to verify our results and Kepler data 
shall be notably useful for such a study. Overall, our procedure provides a key tool for exploring the CoRoT and Kepler databases to 
identify and characterize stellar cycle variability. 

Key words. Physical data and processes: dynamo - methods: data analysis - techniques: photometric - Sun: activity - Sun: magnetic 
fields - Sun: rotation 


1. Introduction 


Dynamo action occurs in the inner stellar layers, but its ef¬ 
fects are clearly observed in and above the photosphere, mostly 
thr ough the magn etic a ctivity cycles. Since the pion eering works 
by IWilsonI (Il978ll and iBaliunas & Vaugh'^ (Il985l) . which pro¬ 
duced the first large sample of magnetic cycles and rota¬ 
tion information validating stellar dynamo models, an increas¬ 
ing number of studies have been dedicated to determining 
of empirical relations between the rotation period Prot, spec¬ 
tral type, cycle lengths, and other stellar parameters in an at- 
tempt to un d erstan d the mechanisms of ste ll ar dynamos (e.g., 
Noves et ^ 1984); Saar & Baliunasl 1992: Baliunas & Sooi 


1995t Brandenburg et alj 11998: Saar & Br andenburg 199* 


a 


Lorente & MontesinosI 20051 iBohm-Vitensel 2007 ). From these 
studies we know now that late-type stars exhibit several chro¬ 
mospheric variability behaviors, ranging from regular variations 
with multi-year periods, similar to the solar cycles, to irregu¬ 
lar cycles with no clear pattern and to stars with no long-term 
variability (IBaliunas & SoonI 119951) . In addition, other studies 
found shown that short stellar magnetic activity cycles, rang- 


magnetic actmty cycles, rang¬ 
ing from 1 to 3 years (iFares et alil2009l: iMetcalfe et al.l ItoTG : 
ISanz-Forcada et ani2013l) . Many authors haye focused on the 
magnetic cycle properties of the Sun; based on sunspot behay- 


* The CoRoT space mission was developed and is operated by the 
French space agency CNES, with the participation of ESA’s RSSD and 
Science Programmes, Austria, Belgium, Brazil, Germany, and Spain. 


ior (see iHathawavl I toToI for a comprehensive review). Gener¬ 
ally, the temporal behavior of sunspots is analyzed using Fo urier 
techniques (e.g.. lMacDonaldlll98^ [Kane & Trivedilll991h . but 
some investigators have employed analytical methods that are 
not based on com binations of periodic sine functions (e.g., 
Mundte^^ 19911) and som e have used wavelet analysis (e.g.. 


Willson & MordvinovI 19991) . 


The amplitude variation of stellar photometric data is still 
poorly studied mainly due to observation length and photomet¬ 
ric sensitivity limitations of the observations. Thanks to the Ke¬ 
pler and CoRoT databases, we now have a unique set of data 
with which to study the nature of these variations. Recent stud¬ 
ies have shown that the photometric variabilities observed by 
Kepler and in data from the Sun behave similarly on all time 
scales, in the sense that the variability indices {Rmr, defined as 
the photometric intensity difference between the 5% and 95% 
value s that contain different v ariability information) are compa¬ 
rable (lBasrietal.ll201lL 12011 . 

More recently. iMathur et al.l (1201 4l) have investigated/Tep/er 
light curves (FCs) of F-type stars searching for indicators of 
magnetic activity cycles, performing a time-frequency analysis 
and using the wavelet procedure. In similar fashion, the standard 
deviation taking in b oxes along of FC has been reported as a 
good activity proxy (iGarcfa et al.l 1201 41) . As revealed by these 
studies, the short and long term behavior of the stellar photo¬ 
metric amplitude variation becomes increasingly important and 
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Fig. 1. The Total Solar Irradiance LC (top-left panel), the prewhitened TSILC (middle-left panels), the amplitude variation A(t) obtained from the 
prewhitened TSI LC (bottom-left panel), and Lomb-Scargle periodograms co rresponding to eac h LC (right panels). In the left panels, the dashed 
vertical lines mark the beginning of each solar cycle, based on review from iHathawa^ l l2010h . In the bottom-left panel, the A{t) time series is 
depicted for Nbcx = 2 (black solid curve) and Ni,„x = 5 (red solid curve), whereas in the bottom-right panel the Lomb-Scargle periodogram of A{1) 
is shown for each Nt,„x value with their respective colors. 


the actual databases allow us to study this relevant attribute of a 
variety of star families. 

This work aims to present an approach to identify and char¬ 
acterize stellar cycles based on LC variations from data collected 
by current photometric space missions. The procedure is then 
applied to a sample of main-sequence CoRoT stars. Section 2 
presents a brief description of the method used to determine 
these cycles, including an application to the Sun’s Total Solar Ir¬ 
radiance. The results for CoRoT stars are presented in Sect. 13.21 
with conclusions in Sect. H] 


2. A search for stellar variability cycles 

Solar cycles have been studi ed from the flux level of solar LCs 
such as VIRGO data (e.g.. lHathawavll20IOl) . Activity-related 
cycles have typically been measured in other stars from the 
photometric variation of activity indexes related with Ca II 
emission. Several works used photometric data from the Mt 


19781 Baliunas & Soonll995USoon et al.ll999l 

Olah et al.l2000l 

Messina & Guinanll2002l; iBerdvugina & Henrv 

l2007HHall et alJ 


20071) . Now CoRoT and Kepler data provide relatively long-term 


LCs with the best photometric sensitivity and temporal resolu¬ 
tion ever. Although long-term variabilities in CoRoT and Kepler 
LCs are affected by instrumental long-term trends that hinder 
physical information, variations such as amplitude are properly 
ret ained in the d i fferen tial LCs. 

iGarcfa et al.l (1201 Ol) presented for the first time a method to 
study long-term stellar cycles using the CoRoT LC of HD 49933 
provided in the astero-field channel with a time sampling of 32 s. 
Based on asteroseismology studies, these authors detected long¬ 
term variations of the amplitude, frequency shift, and standard 
deviation of the modulation produced by the acoustic-mode en¬ 
velope. These variations were obtained from subseries of 30 
days, shifted every 15 day s . A si milar procedure was recently 
performed bv iMathur et al.l (1201 3h to study the stars HD 49385, 
HD 181420, and HD 52265 using CoRoT LCs combined with 
NARVAL (e.g.. I Aurierel2003T) spectropolarimetry data. More re¬ 


cently, iMathuTeril] (I 2 OI 4 I) detected long-term variations in Ke¬ 
pler LCs of 22 F-stars by considering time-series of different ac¬ 
tivity proxies, such as the flux standard deviation and temporal 
sections of the wavelet power spectrurrQ. In particular, the stan¬ 
dard deviation of LC subseries has been suggested to be a good 
activity proxy, named Sph^k, where k is a multiplier of Prot to 
define the leng th of the subseries u sed to measure the magnetic 
proxy. As such. lGarcfa et al.l(l2014l) also considered this proxy in 
a study o f activity cycles dete cted in 540 Kepler puls ating sun¬ 
like stars. iMathuret al.l(l2014l) and iGarcfa et al.l(l2014l) have sug¬ 
gested that computing S ph^k within time boxes of 5 y. Prot, namely 
Sph,k= 5 , provides a proper detection of Sun-like activity cycles. 
The activity proxies of those studies are related with the overall 
LC variability, which for those Kepler sources are dominated by 
the rotational modulation. 

In the present work, we consider an appro ach somewhat sim¬ 
ilar to that performed bv iMathur et al.l (l2014l) . However, we sug¬ 
gest as an appropriate activity proxy the amplitude of the rota¬ 
tional modulation computed within random time boxes, as de¬ 
scribed below. This proxy was used to compute stellar cycles for 
a sample of 16 CoRoT stars, with a first application to the Total 
Solar Irradiance (TSI) LC. The proxy is clearly valid because, 
for instance, the variability amplitude is somewhat related with 
the LC standard deviation, as well as with the wavelet power 
spectrum. In particular, the wavelet analysis has the limitation 
of requiring a continued time covering to be properly applied. 
In contrast, our method is an important application for LCs with 
long-gaps, as detailed in Sect 3.2. 

To further elaborate, our procedure is described by the fol¬ 
lowing steps: 

1. For homogeneity, either a TSI or a CoRoT LC is analyzed in 

normalized flux units, obtained by dividing the original LC 

by its entire flux average. 

* This spectrum is somewhat comparable to a time-frequency analy¬ 
sis and consists of a correlation between the time series and a mother 
wave l et given with diffe rent periods along time dTorrence & Comp^ 
ll998l:lMathur et al.l2014l) . 
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Table 1. CoRoT stars analyzed in this work. 


CoRoT parameters 

Our parameters 

CoRoT ID 

ST 

LC 

B 

(mag) 

V 

(mag) 

id 

Pmt 

(days) 

Pcyc Sequence" 

(days) 

Q 

102712791 

G4 

IV 

13.943 

13.666 

a 

0.96 + 0.03 

33 + 3 

IS 

~ 1 

102777006 

G4 

V 

16.542 

15.365 

b 

1.33 + 0.02 

101 + 8 

I 

~ 1 

102750723 

G1 

V 

15.162 

14.019 

c 

1.44 + 0.02 

106 + 7 

I 

~ 1 

102721955 

KO 

IV 

14.641 

14.210 

d 

2.17 + 0.06 

187 + 20 

I 

0.85 

102752622 

G2 

V 

14.682 

14.208 

e 

2.33 + 0.07 

346 + 24 

AI 

0.59 

102780281 

F8 

V 

15.889 

14.585 

f 

3.0+ 0.1 

110 + 8 

IS 

~ 1 

102726103 

G4 

IV 

15.096 

14.329 

g 

3.7+ 0.1 

277 + 21 

I 

0.64 

102770332 

G7 

V 

17.644 

15.589 

h 

4.2+ 0.1 

424 + 41 

I 

0.42 

102770893 

G5 

V 

15.169 

14.295 


4.3 + 0.2 

427 + 45 

I 

0.42 

102749950 

G1 

V 

14.883 

14.226 

j 

5.4+ 0.2 

408 + 26 

I 

0.42 

102758108 

G5 

V 

15.098 

14.457 

k 

6.1 +0.2 

614 + 55 

I 

0.34 

102754736 

K2 

V 

14.825 

14.345 

1 

6.9+ 0.3 

117 + 8 

S 

~ 1 

102723038 

G8 

V 

15.881 

14.477 

m 

8.6 + 0.5 

210 + 7 

S 

0.78 

102720703 

K3 

V 

14.376 

13.862 

n 

10.2 + 0.6 

650+ 130 

I 

0.34 

102778595 

KO 

IV 

15.857 

14.700 

0 

11.8 + 0.7 

201 + 15 

S 

0.82 

102738457 

G1 

V 

15.639 

15.047 

p 

12.9 + 0.6 

239 + 22 

S 

0.78 


Note. 

“ Sequence refers to the likely classification of activity level in the versus P^oi diagram (see Sect. 13.4b based on previous studies (e.g., 
lBdhm-Vitensell20(y^ : ISaar & BrandenburglfT9^ : lBrandenburg et alJIlQQ^l (A of active and I of inactive) and in our new results (S of short cycle 
periods). AI or IS means an intermediary location between A and I or I and S. 


2. The LC is then prewhitened. For CoRoT LCs, multiple runs 
are then combined to provide full time series for each tar¬ 
get. This prewhitening must be performed on CoRoT LCs 
because these data are subject to long-term trends that mask 
physical signatures at long periods. As such, the TSI data 
are also prewhitened to simulate the CoRoT LC limitations. 
Details regarding this prewhitening will be discussed below. 

3. We compute Proi as being the main peak of a Lomb- 
Scargle periodogram (iLombll 19761: IScar ^ Il982h of the en¬ 
tire prewhitened LC. The periodogram frequency range is 
dehned between /o = 2/7'tot, where Ttot is the total baseline 
of the observations in days, and - 10 d“', with 10^ el¬ 
ements. We r ehned Pmt by ana lyzing its dispersion dehned 
by Eq. (2) of iDworetskvl (Il983h as a function of period. We 
computed a Gaussian ht with this function in such way that 
its minimum provided the rehned Prot value and its full width 
at half maximum (FWHM) was assumed to be the Prot uncer¬ 
tainty. More details about this Pr„t calculation are discussed 
below. 

4. The LC is analyzed within random time boxes with sizes 
given by NboxProt with initial times f,-, where Nbox is the num¬ 
ber of rotation cycles within each box. The choice of the Nbox 
value is discussed below. 

5. For each box, the phase diagram is computed for a pe- 
riod equal to Pmt- Then, a harmonic ht is computed as in 
iDe Medeiros et all (1201 3h . given by the following: 

M 

y(f) = ^ [a j sin (27r/i jt) -H h; cos (27r/i ;f)], (1) 

1=1 

where /i = 1 IProt, aj and hj are the Fourier coefhcients, M 
is the number of harmonics, and t is the time. For this work, 
we set M - 2, based on a compromise between a reason¬ 
able ht and computation time. From this ht, the amplitude 
A, is computed as the difference between the maximum and 
minimum values of y{t). 


6 . The above steps are repeated in N iterations to obtain a suit¬ 
able prohle of the amplitude variation over time. For this 
work, we assign N = lOTtotIProt, which corresponds to ap¬ 
proximately 10 amplitude measurements per rotation cycle. 
This number is chosen in a compromise between the stability 
of the results and computation time. 

7. Next, a time series A{t) is dehned as being the set of A,- values 
corresponding to each f, value. 

8 . Finally, the long-term cycle period {Pcyc) is computed as the 
main peak of a Lomb-Scargle periodogram of the A(f) time 
series. For this periodogram, the frequency range is dehned 
between/o = 2/7’,o, and/;v = 11 Prot, with 10^ elements. We 
also rehned and estimated its uncertainty by following 
the same procedure as described for Pr„t in step 3. 

Regarding the second step, the prewhitening was tested in 
various ways and provided stable results. Overall, for CoRoT 
data, the prewhitening was obtained by dividing the LC of 
each run by a third-order polynomial ht, as typically per- 
formed in the litera ture (e .g., lBasri et alJl201 ItlAffer et alJJoiTl : 
iDe Medeiros et al]l2013h . For TSI data, the prewhitening was 
obtained by diving each LC by a boxcar smoothed version of 
itself. The time box in this case can be dehned within approx¬ 
imately 100-200 days, corresponding to the typical time win¬ 
dows of the CoRoT runs (of ~150 days). Figure [1] (middle- 
left panel) illustrates a prewhitened TSI LC for a time box 
of 150 days. 

As to the third step, those simple period estimations are 
valid for approximately 95% of all ta rgets with LCs showing ro - 
tational modulation, as estimated in iDe Medeiros et"^ (l2013h . 
The remaining 5% typically comprise misdetections of aliases 
for stars with main spots at opposite longitudes, as explained in 
the referred work. We also avoided the surface granulation sig¬ 
nal, which in the TSI LC is a noisy signature distributed at sev¬ 
eral frequencies of the power spectrum (see, e.g., Agrain et al. 
2004 and references therein). This so-called Solar background 
is mostly concentrated at low frequencies (typically <0.1 d *) 
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Fig. 2. CoRoT LCs and the A(t) time series for typical sources in our final sample. For each star we show the combined LC (upper panel), with the 
corresponding A(t) time series (bottom panel). The ID and computed and are given above each panel. 


where the Solar Prot lies. However, that background does not 
hinder the Proi determination because it has periodogram pow¬ 
ers considerably weaker than the p eaks produced by the rotating 
spots. Although we can infer from I Ai grain et akl (l2004l) that the 
Solar background may be related with the activity long-term cy¬ 
cles, exploring such a relation is not the scope of our work. As 
such, for the CoRoT LCs, we assumed the granulation signal 
to be a noise, which was avoided b y selecting LCs w i th hig h 
amplitude-to-noise levels, based on lOe Medeiros et al] (1201 3b : 
see also Sect. 3.2. 


In the fourth step, the choice for random f, values is not 
mandatory. However, using such an irregular box sampling min¬ 
imizes possible window function problems that could be pro¬ 
duced by regular bins. In addition, the Nbox value is not unique. 
We used Nhox = 2 to have at least two cycles in each box for the 
calculation of the phase diagram and its amplitude and, at the 
same time, to have a high sampling in A{t). For validation, we 
verified that varying Nbox with slightly greater values provides 
similar results. Figure[T](bottom panels) illustrates the main data 
in our analyses for Nbox = 2 (black solid curves) and 5 (red 
solid curves). Although the amplitude variation plotted in the 
bottom-left panel is smoother for Nbox = 5 than for Nbox - 2, 
both values produce similar periodogram peaks, as shown in the 
bottom-right panels. Therefore, the final results (see Sect. [3j for 
long time variations are stable over different Nbox values. 


3. Results 


3.1. Revisiting the Totai Soiar Irradiance tight curve 


TSI data h a ve bee n collected from space over the past 35 years. 
iHathawa^ (l2010ll presents a chronological list of the instru¬ 
ments that have performed these measurements, with a review 
of the main results achieved using these data, based on compos- 
ite TSI time se r ies constructed using d i fferen t procedures (e.g., 
WillsonI I1997I: IWillson & MordvinovI I1999I: iFrohlich & Leanl 


19981). Here, we revisit the TSI LC in normalized flux units to 


validate the application of our new approach to CoRoT LCs. 

Figure[T]shows the original and prewhitened TSI LCs (upper 
and middle left panels, respectively). The Lomb-Scargle peri¬ 
odogram of the original LC (see top-right panel) shows a promi¬ 
nent peak, which is the well-known 11-yr solar cycle. However, 
as expected, this period disappears in the periodogram of the 
prewhitened LC (see middle-right panel). Therefore, long period 
could not be properly detected in other stars from the CoRoT 
flux variations because of the aforementioned long-term trends. 
In our new approach, the amplitude variation (bottom-left panel) 
was obtained from the prewhitened LC and, in this case, the 11- 
yr solar cycle was confidently recovered in the power spectrum 
(see bottom right-panel). Therefore, our method is suitable to 
search for long-term variability cycles with long time span data. 


In Fig. [T] the power spectrum of the original TSI LC (see 
top-right panel) is dominated by long-term variations (including 
11 -yr solar cycle), whereas the periodogram of the prewhitened 
LC is dominated by short-term variability (see middle-right 
panel). As such, we first computed Pyot from the prewhitened 
TSI LC (middle-left panel), as described in Sect. |2] The Lomb- 
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ID - 102712791 = 0.09(yr) 
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Fig. 3. Phase diagrams from A{t) for the stars referred to in Fig.|2] with 
the corresponding ID indicated in the header. 


Scargle periodogram (middle-right panel) yielded a main pe¬ 
riod of 30.6 days, which cor respond to the sola r surface Prot at 
its active latitudes (see, e.g.. iLanza et al.ll2003h . From the A{t) 
time series (lower-left panel), the power spectrum (lower-right 
panel) recovers well-known long-term periods from the litera¬ 
ture: 11.6 + 0.9 yr, 8.3 + 0.7 yr, 1.2 + 0.4 yr, 2.4 + 0.6 yr. 

The periods of 1.2 yr and 2.4 yr computed here are compati¬ 
ble with the typical time scales o f 1.1 yr and 2.8 yr , respectively , 
of the so-called Gnevyshev gaps (iGne wshevi 1 963l 1 1 967i (197% . 
These gaps, which are relatively short-term decreases in the solar 
activity, can be seen in the A{t) time series (Fig. [T] bottom-left 
panel), especially considering the red solid curve. For instance, 
during the solar maximum phase of cycle 22, there is a promi¬ 
nent short-term decrease throughout the years 1989-1991 and 
several fainter decreases are also noticeable during other maxi¬ 
mum phases. The 1.2 yr period computed here is also compatible 
with a quasi-periodicity of approximately 1.3 yr, which is related 
to the stochastic processes of active region emergences produced 
by the solar magnetic d ipole moment and open magnetic flux 
(IWang & Sheelevll2003h . Another periodicity of approximately 
1.0 yr is also produced by dynamical processes of differential 
rotati on, meridiona l flow, and sup e rgranu le diffusion (iHathawavl 
I 2 OIOI) . In contrast, iFletcher et al.l (1201 Ol) suggests that a period 
of 2.0 yr, which is most prominent during the solar maximum 
phases, may be related to the acoustic modes of the main 11-yr 
solar cycle. Hence, several phenomena are related to these dif¬ 
ferent activity cycles identified in the TSI LC, such as the main 
11-yr cycle, as well as the shorter cycles of 1.2 + 0.4 yr and 
2.4 + 0.6 yr computed here (1-2-yr cycles henceforth). There¬ 
fore, detecting either 11-yr- or 1-2-yr-like cycles in other stars 
is an important goal for studying activity cycles. 

3.2. Stellar cycles from CoRoT light curves 

To identify stellar cycles in other stars, we analyzed 877 CoRoT 
LCs from the exoplanet channel, which are publicly available 
in a reduced form at obtained by the so- called CoRoT N1-N2 
pipelineQ (see, e.g.. I Auvergne et ^l2009h . We selected the LCs 
that provide the longest combined time coverage, namely those 

^ http://idoc-corot.ias.u-psud.fr/ 
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Fig. 4. Typical examples of the simulated LCs produced from the TSI 
data. From the top to bottom panels, these LCs have time coverages of 
approximately 800, 4000, and 8000 days. 


comprising the observing runs IRaOl, LRaOl, and LRa06, which 
have typical time spans of 57, 131 and 76 days, respectively. 
Thus, these LCs combined for each star have coverages of ~ 
264 days and a total time span of ~ 1880 days (~ 5.1 years). 
From visual inspection we identified 16 FGK type stars with 
clear rot ation signatures and an am plitude-to-noise ratio, as com¬ 
puted in iDe Medeiros et^ (1201 3h . greater than 2. These 16 se¬ 
lected stars are listed in Table 1. In addition, a cross-check with 
the literature showed that 14 stars of this sub-sample were al¬ 
ready known as p hotometric variables (iDe Medeiros et all2013l : 
lAffer et £012012h . Figure illustrates typical LCs of our sub¬ 
sample (upper panels) with their respective A(t) time series (bot¬ 
tom panels). The computed Pcyc and Prot values of the entire 
sample of 16 stars are listed in Table 1, with Prot and Pcyc ranging 
from 0.956 to 12.903 days and from 0.09 to 1.77 years, respec¬ 
tively. The table also provides the recovery fractions Q obtained 
from simulations, which are described in Sect. 13.31 as well as 
target identifiers (id) and their evolutionary sequences A, I or S, 
which are considered in Sect. 13.41 

Numerically, the A(t) time series were analyzed in this work 
using a similar technique as that considered in the analysis of ra¬ 
dial velocity time series RV(f). Although the physical informa¬ 
tion of A{f) is obviously different than that provided by RV{t), 
the numerical techniques involved also consist of identifying a 
persistent signal in the phase diagram. Therefore, we only need 
data points in different phases and levels of activity to deter¬ 
mine the stellar cycle periods from A(f) (see Sec. 13.31) . Figure[3] 
shows the phase diagrams obtained from the four LCs displayed 
in Fig.|2] These phase diagrams exhibit a well-defined signature 
by considering data points that can be distributed along sepa¬ 
rated temporal windows (i.e., time series with gaps). All periods 
w ere computed w ith a false alarm probabi l ity (FA P), as defined 
in lScargle](ll982l) and in iHorne & Baliuna'sl(ll986l) . less than 0.01 
(i.e., a significance level > 99%). We verified from simulations 
(see Sect. 13.31) that the periods compatible with long-term cycles 
have a mean \og{FAP) ~ -10, whereas incompatible periods 

have a mean log(F’AP)-6. Despite this difference, in any case 

the FAP values are lower than 0.01. Hence, these low FAP values 
indicate that the periods detected in the periodogram maxima are 
statistically acceptable; however, they can have different natures 
than actual cycles. We present below an approach to verify the 
confidence of these periods regarding their compatibility with 
long-term cycles and considering the CoRoT data limitations. 
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Fig. 5. Recovery fraction Q as a function of TcoverIPcyc obtained from 
10'* simulations where random long-gaps were inserted in the TSI LCs 
(see text for details). 

3.3. On the influence of gaps 

The entire LC of each target being combined from multiple runs 
are not continuous for each target, as described in Sect. |2 Thus, 
the full LCs may exhibit long gaps that may hinder a proper de¬ 
tection of solar-like cycles such as the 11-yr or the 1-2-yr cycles 
described in Sect. 13.11 To test how these gaps may affect the de¬ 
tection of variability cycles in the CoRoT LCs, we performed 
several simulations by inserting synthetic gaps into the TSI LCs. 

Specifically, 10^ LCs were automatically treated and ana¬ 
lyzed as follows. In each simulated LC, which originally cover 
35 yr, a certain number of random gaps with durations ranging 
from 200 to 4000 days was inserted such that time coverages of 
400 to 8000 days remained in the simulated LC. Figure |4] shows 
typical examples of the resulting simulations with different time 
coverages. From the entire set of simulations, we computed the 
recovery fraction Q of the 11-yr cycle as a function of the time 
coverage to Pcyc ratio, or TcoverIPcyc, as depicted by the black 
histogram in Fig. |5] This calculation was performed within Pcyc 
bins of 0.25 yr by considering that the 11-yr cycle was properly 
recovered if the main periodogram peak of the simulated LC fell 
within 11.6 + 0.9 yr, the central value and uncertainty computed 
from the original TSI LC. A low value of Q does not necessarily 
means a misdetection, but it means a less reliable measurement 
of the period. Indeed, the recovery fraction increases if a higher 
uncertainty is assumed for the 11-yr cycle period, as is the case 
of some cycles detected in the sample of table [T] 

Based on Fig.|5]and on the re lation between Pcyc and Pr„t ob¬ 
tained by iBohm-Vitensel (l2007l) . we expect to detect 11-yr-like 
cycles with high confidence (Q > 75%) for Proi 5 2. This covers 
the first three or five targets of table[T] Possible 11-yr-like cycles 
were detected in some targets with Prot > 2 d. Even though these 
cycles have low Q values, we consider them for a discussion in 
Sect. l3.4l about their overall distribution. The simulations predict 
high Q values for 1-2-yr-like cycles, suggesting a high confi¬ 
dence for those cycles. However, one should be aware that these 
cases were assumed as being short versions of the 11-yr cycle 
within gapped LCs. Estimating Q values for CoRoT data based 
the own 1-2-yr solar cycles is not obvious because the detection 
of these cycles in the TSI data is hindered by the strong superim¬ 
posed 11-yr cycle. In addition, the 1-2-yr cycles have low signals 
in the TSI LC, in contrast with the possible 1-2-yr-like cycles de¬ 
tected in the CoRoT LCs. Because of these differences between 
the TSI and CoRoT signals a high Q value can only tell us that 



0,0 0.5 1,0 1.5 2.0 

log Prot(dQys) 


Fig. 6. The P^yc (in years) as a function of the Pro, (in days) on loga¬ 
rithmic scales, with our data represented by red (Q > 0 .75) and blue 
(Q < 0.75) circles, da ta from ISaar & Brandenburgl7l999h by black tri¬ 
angles and data from iLovis et all l l201 iF by grey circles. The squares 
indicate the sources with secondary periods for s ome stars on the A 
sequence according to I Saar & Brandenburd lll999h . The Sun magnetic 
cycle periods of 11.6 yr (black Sun symbol) and 1.1 yr (gray Sun sym¬ 
bol) are also represented. The Pcyc vers us Prot empiric a l relat ions (full 
lines A and I sequences) proposed by iBbhm-Vitensd ( 120071) are also 
displayed. The empirical relations for shorter term stellar magnetic cy¬ 
cle periods (dashed line S sequence) were obtained by shifting the I 
sequences by a factor of 1/4. 


a cycle was well-measured because of their long time-coverage, 
without taking into account other limitations such as the possible 
hindrance produced by a strong superimposed signal. 

3.4. On the physical nature of the cycle periods 

Previous works (e.g..lBohm-Vitensj|2007l : ISaar & Brandenburg 

1 19991: iBrandenburg et al.lll998h ha^ identified that Pcyc and Prot 

(as well as the Rossby number) are correlated one another show¬ 
ing an evolutionary behavior with two main sequences: A (ac¬ 
tive) and I (inactive), which are compose d of young and o l d stars , 
respectively. Following the approach of iBohm-Vitensd (l2007h . 
we present in Fig. |6]the relation between Prot and Pcyc in loga¬ 
rithmic scale for our CoRoT sample. These are plotted with data 
from previous works for co mparison, as well as w ith the known 
parallel A and I sequences. iBohm-Vitensel (l2007h also detected 
a split into the Aa and Ab sequences, which are differed from 
the I sequence by a factor of ~6 and ~4 (for a given Prot), re¬ 
spectively. We consider in Fig. |6] the sequence Ab. We suggest 
here to include a third sequence, namely the S sequ ence, based 
on the data from CoRoT combined with those from iLovis et ^ 
dMni). This S sequence, depicted in Fig.|6] parallels the A and I 
sequences and it differs from the I sequence by a factor of ~ 1/4 
(for a given Prot). 

The compatibility of the CoR oT data dist r ibutio n as an ex¬ 
tension of the distribution from ILovis et ^ (1201 ih in Fig. |6] 
supports a validation of the S sequence. However, more data 
is needed for such a final confirmation. Besides, establishing a 
physical interpretation for this sequence is difficult in the present 
study because of the limited number of targets and of the CoRoT 
data limitations. Longer and higher-quality time series provided 
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by the Kepler mission a s already done for a sub -sample of stars 
by Mathur et alJ (l2014l) . IVida et all (1201 4h and lArkhvpov et alJ 
(1201 5I) shall help greatly a more in-depth study of the S sequence 
. Currently we suggest that the S sequence is possibly physical, 
but it can come from a bias produced by the CoRoT time window 
limitation, thus exhibiting 1-2-yr-like cycles of stars that would 
actually lie in the A or I sequence. As such, these short cycles 
could be physically related with the Solar Gnevyshev gaps. As a 
second possib ility the S sequence co uld fulfill a third evolution¬ 
ary region, as iBohm-Vitensel (l2007h suggested for the A and I 
sequences. These possibilities shall be explored in further stud- 


There are some stars lying in intermediary regions between 
two sequences, as depicted in table [1] In particular, targets a and 
/ lie between the I and S sequences and have high Q values. 
This suggests that their locations in Fig. 0 are confident and they 
could be experiencing a transition between two sequences, as 
we could say about the Sun, which lies between the A and I 
sequences. A similar case to the Sun could be target e. However, 
we cannot consider its location with confidence because of its 
relatively low Q value. 


4. Conclusions and Future Works 

In this study we have shown that the A(f) time series, which ex¬ 
presses the behavior of the amplitude variation as a function of 
time, seems to enable an important procedure to compute stellar 
cycles from photometric LCs. In this context, to check our pro¬ 
cedure, we have shown that the variability periods for the Sun 
from the A(f) time series agree with those computed in previous 
works; that used different approaches. In addition, the shape of 
the A{t) time series allows us to determine the beginning of the 
grand minimum and examine the global behavior of the solar 
magnetic cycle. 

From the A{t) time series, we have computed variability cy¬ 
cles and rotation periods for 16 FGK CoRoT stars. The behavior 
of the obtained cycle periods in the Pcyc versus distribution 
follows a similar relation found in other studies for other families 
of stars, which suggests a physical nature of the identified stellar 
cycles. Simulations from Solar TSI data indicate that about a half 
of our sample has a recovery fraction Q > 0.75 and, therefore, 
further study is needed for more conclusive results. The confir¬ 
mation of the Pcyc versus Prot empirical relation identified in this 
work combined with those identified in previous studies (e.g., 
lBohm-Vitensell2()(y^ may allow us to develop a more robust ex¬ 
planation for their physical nature. 

To date, the long-term stellar variability of stellar LCs, such 
as those collected by the CoRoT and Kepler space missions, re¬ 
mains largely unexplored. In many cases these stellar LCs may 
contain large gaps or instrumental trends that hinder identifying 
long-term physical variations. The current procedure, namely the 
A{t) time series approach, offers a clear methodology to identify 
and compute variability cycles from the photometric database of 
these missions. As suc h, the longer time s pan provided by Kepler 
and PLATO missions (iRauer et al.ll2014l) . shall allow us to com¬ 
pute more reliable Pcyc measurements and thus a more in-depth 
study about stellar cycles. 
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